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Abstract 

A progressive failure analysis methodology has been 
developed for predicting the nonlinear response and 
failure of laminated composite structures. The 
progressive failure analysis uses C 1 plate and shell 
elements based on classical lamination theory to 
calculate the in-plane stresses. Several failure criteria, 
including the maximum strain criterion, Hashin’s 
criterion, and Christensen’s criterion, are used to predict 
the failure mechanisms. The progressive failure 
analysis model is implemented into a general purpose 
finite element code and can predict the damage and 
response of laminated composite structures from initial 
loading to final failure. 

Introdtiction 

Composite materials have been increasingly used in 
aerospace and automotive applications over the last 
three decades and have seen a dramatic increase in usage 
in non-aerospace products in the past five years. The 
use of composite materials is very attractive because of 
their outstanding strength, stiffness, and light-weight 
properties. An additional advantage of using 
composites is the ability to tailor the stiffness and 
strength to specific design loads. Since most composite 
materials exhibit brittle failure with little or no 
ductility, as offered by metals, the behavior of the 
composite structure must be understood, and analysis 
methods to predict and propagate the failure need to be 
developed. Laminated composite structures can develop 
local failures or exhibit local damage such as matrix 
cracks, fiber breakage, fiber-matrix debonds, and 
delaminations under normal operating conditions which 
may contribute to their premature failure. The ability 
to predict the initiation and growth of such damage is 
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essential for predicting the performance of composite 
structures and developing reliable, safe designs which 
exploit the advantages offered by composite materials. 
Hence, the need for a reliable methodology for 
predicting failure initiation and propagation in laminated 
composite structures is of great importance. 

In recent years, the progression of damage in 
composite laminates has been a focus of extensive 
research. Ochoa and Reddy present an excellent 
discussion of progressive failure. 1 Reddy and Pandey 2 
developed a finite element procedure based on first-order 
shear deformation theory for first-ply failure analysis of 
laminated composite plates subjected to in-plane and/or 
transverse loads. Pandey and Reddy 3 extended their 
earlier work on first-ply failure of two-dimensional 
laminated composites to include a progressive failure 
analysis capability. However, only a linear finite 
element analysis was performed. Reddy and Reddy 4 
calculated and compared the first-ply failure loads 
obtained by using both linear and nonlinear finite 
element analyses on composite plates. The differences 
between the linear and nonlinear failure loads were found 
to be large for the cases of transverse loading, and were 
considerably less for the cases of in-plane (tensile) 
loading. Ochoa and Engblom 5 presented a progressive 
failure analysis for composite laminates in uniaxial 
tension using a higher-order plate theory with shear 
deformable elements. Analyses were performed on a 
plate subjected to uniaxial tension and to four-point 
bending. However, comparisons with experimental 
results were not provided. Engelstad, Reddy, and 
Knight 6 investigated the postbuckling response and 
failure prediction of flat composite unstiffened panels 
loaded in axial compression using 9-node shear 
deformable elements. Good correlation between the 
experimentally obtained and analytically predicted 
postbuckling responses was observed for deflections and 
surface strains. Only the maximum stress criterion and 
the Tsai-Wu criterion were used in the failure 
predictions. Chamis et al. 7 presented a damage tolerance 
model for pressurized cylinders which included damage 
detection, progression, and failure. Failure criteria and 
delamination models were assessed and material 
properties reconstituted. Only analytical simulations 
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were presented. Hwang and Sun 8 performed a failure 
analysis of laminated composites by using an iterative 
three-dimensional finite element method. The three- 
dimensional analytical results agree favorably with the 
experimental results for notched and unnotched 
specimens loaded in tension. However, the analytical 
predictions underestimated the experimental results for 
angle-plied laminates with holes. Chang and Chang 9 
developed a progressive failure damage model for 
laminated composites containing stress concentrations. 
They applied this progressive failure method to bolted 
composite joints and to a laminated composite plate 
containing a hole. Comparisons were made to 
experimental results in these studies and reasonable 
correlation to the data was reported. In summary, 
whereas much research has been performed on the 
progressive failure of structural components undergoing 
linear deformations, limited research including nonlinear 
deformations has been done. 

The objective of this paper is to examine the use of 
progressive failure analysis in the nonlinear deformation 
regime. A progressive failure methodology is developed 
and implemented into a general purpose finite element 
code and validated by comparing analytical predictions 
of progressive failure of geometrically nonlinearly 
deformed composite structures with experimental data 
Different formulation methods for detecting failure am 
also compared and assessed. This effort incorporates the 
failure detection and material degradation models of 
different researchers into a single computational 
structural mechanics framework. As such, different 
approaches can be readily compared and new failure and 
materials models incorporated. 

Progressive Failure Methodology 

A typical methodology for a progressive failure 
analysis is illustrated in Figure 1. At each load step, a 
nonlinear analysis is performed until a converged 
solution is obtained assuming no changes in the 
material model. Then using this equilibrium state, the 
stresses within each lamina are determined from the 
nonlinear analysis solution at this step. These stresses 
are then compared with material allowables or used to 
assess certain failure theories. A survey of failure 
criteria is presented by Tan. 10 If a lamina failure is 
detected, as indicated by a failure criterion, the lamina 
properties are changed according to a particular 
degradation method. Then, the initial nonlinear 
solution no longer corresponds to an equilibrium state 
and equilibrium of the structure needs to be re- 
established utilizing the modified lamina properties for 


the failed lamina while maintaining the current load 
level. This iterative process of obtaining nonlinear 
equilibrium solutions each time a local material model 
is changed is continued until no additional lamina 
failures are detected. The load step is then incremented 
until catastrophic failure of the structure is detected. 

Therefore, typical progressive failure analysis 
methods involve five key features. First, a nonlinear 
analysis capability is needed in order to establish 
equilibrium. Second, an accurate stress recovery 

procedure is needed in order to establish the local lamina 
stress state. Third, failure detection criteria are needed in 
order to detect or sense local lamina failure which may 
also indicate the nature of the failure. Fourth, material 
degradation or damage models are needed in order to 
propagate the failure and establish new estimates for the 
local material properties. Finally, a procedure to re- 
establish equilibrium after modifying local lamina 
properties is needed. This research focuses on the last 
four features since nonlinear analysis procedures are 
already well established. The progressive failure 
methodology developed in this paper is an extension of 
work performed by Pifko 11,12 of Northrop-Gnmiman 
Aerospace Corporation.. 

To implement and perform the progressive failure 
methodology described earlier, a structural analysis 
software system is needed. The framework to 

accomplish this task is a structural analysis software 
system called COMET (Computational MEchanics 
Testbed) 13,14 developed jointly between NASA Langley 
Research Center and Lockheed Palo Alto Research 
Laboratory. COMET’S capabilities include linear and 
nonlinear stress analyses of large-scale built-up 
structures, transient dynamic analyses, and eigenvalue 
analyses. 

The progressive failure analysis methodology uses 
C 1 plate and shell elements based on classical 
lamination theory to calculate the in-plane stresses. 
The nonlinear Green-Lagrange strain-displacement 
relations are used in the element formulation, and large 
rotations are treated through the element-independent 
corotational formulation in COMET. The progressive 
failure methodology implemented in COMET 
accommodates the maximum strain criterion, 
Christensen’s criterion, 15 and Hashin’s criterion. 16,17 
When a failure is detected, the progressive failure model 
classifies the mode of failure as fiber failure, matrix 
failure, shear failure, or an interaction of failure modes. 
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In the maximum strain criterion, failure is assumed 
to occur if any of the following conditions are satisfied: 

£, > X £t or |e, | > X g( , fiber failure 
e 2 > Y £t or |£ 2 |>K £c matrix failure (1) 
T e shear failure 

where e ly e 2 ^nd y l2 are the lamina strains in the 
material coordinate system and X £t , X £c , Y £t , 
Y £c , and T £ are the material allowable strains denoted as 

X e - critical tensile strain in fiber direction 
€ r 

X £c = critical compressive strain in fiber direction 

K - critical tensile strain in matrix direction 
e T 

Y £c = critical compressive strain in matrix direction 
T e = critical shear strain 

The absolute value sign on y 12 indicates that the 
sign of the shear strain is assumed to not affect the 
failure criterion. The maximum strain criterion is a 
non-interactive failure theory in strain space. Since the 
maximum strain criterion provides different conditions 
for failure, the mode of failure can be identified as either 
fiber failure, matrix failure, or shear failure. 

Hashin 17 proposed a stress-based failure criterion that 
has the ability to predict the modes of failure. The 
observation of failure in unidirectional fibrous 
composites indicates that there are two primary failure 
modes: a fiber mode in which the composite fails due to 
fiber breakage in tension or fiber buckling in 
compression; and a matrix mode in which matrix 
cracking occurs. Since different failure mechanisms 
occur in tension and compression, Hashin further 
subdivided each failure mode into a tension and 
compression mode. The failure modes are summarized 
for the case of plane stress as follows: 

Tensile Fiber Mode, o x > 0 



Compressive Fiber Mode, <J X < 0 



Tensile Matrix Mode, <r 2 > 0 



Compressive Matrix Mode, a 2 < 0 



where <r l9 <r 2 ,and r 12 are the in-plane lamina stresses 
and X Ti X c , Y Ti Y c , and T are the lamina strengths. 
The subscripts T and C in X and Y refer to the normal 
strengths in tension and compression. 

Christensen 15 introduced a quasi-three-dimensional 
lamination theory which accounted for the out-of-plane 
stress terms. In related work, Christensen then 
proposed a strain-based failure criterion which 
distinguished the modes of failure into either fiber 
failure or fiber/matrix interaction failure. The 

corresponding equations for failure are as follows: 

Fiber Failure: 

X Ec Ze x <X Sr ( 6 ) 

Fiber/Matrix Interaction Failure: 

P £ kk + e ij e ij - a2 (7) 

where /? and a are determined from experimental 
failure data and is the deviatoric strain tensor given 
by 

&ij — £jj — ( 8 ) 
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In Christensen's analyses, the two parameters /? and (X 
in Equation (7) are evaluated to fit failure data for tensile 
and compressive failure with no shear stress. 

When a failure is detected, the progressive failure 
model classifies the mode of failure as fiber failure, 
matrix failure, or shear failure. Two material 
degradation models are implemented including 
instantaneous reduction and gradual reduction of the 
material properties for use with the ply-discount 
approach. The material properties which are degraded 
when a failure is detected depend upon the failure mode. 
The solution procedure performs a geometrically 
nonlinear analysis using a modified Newton-Raphson 
algorithm for either applied force or applied 
displacement problems. The procedure is also capable 
of using a global load-stepping algorithm for advancing 
the nonlinear analysis solution during a static analysis. 
Equilibrium is not presently re-established in the 
nonlinear analysis procedure implemented in COMET. 
Instead, the strategy is to use small load increments in 
the nonlinear analysis procedure which minimizes the 
effect of not re-establishing equilibrium as shown in 
Figure 1. 

Results 

The progressive failure methodology described earlier 
has been successfully implemented in COMET. Two 
examples are considered for the evaluation of the 
methodology. The first example deals with a composite 
laminate under rail-shear loading. The second example 
is a laminated composite panel subjected to an axial 
compressive load. Numerical results obtained using the 
present progressive failure methodology are compared 
with other published results. Additional problems have 
been solved and are reported by Sleight. 18 

These examples are of interest because of the 
available experimental data and their applicability to 
aircraft structures. The first example is a membrane 
problem with the nonlinear behavior due to material 
failure. The second example involves combined 
membrane and bending behavior and combined material 
degradation and large deflections. 

For the models considered, determining the structural 
response involves several steps using COMET. For the 
model with compressive loading, the first step is to 
perform a linear static analysis. Using the linear static 
analysis, the primary equilibrium path is found which 
produces no out-of-plane deflection. The next step is a 
linear stability analysis to find the point at which the 


primary path will bifurcate to a secondary equilibrium 
path. Along the secondary path, the transverse 
deflections will increase. The third step is forming the 
initial geometric imperfection which serves as a trigger 
for the geometric nonlinearities. For the model with 
shear loading, these first three steps are not required 
since buckling does not occur. For both models, the 
next step is to perform an elastic nonlinear analysis to 
understand how the structure behaves without any 
material failures. Finally, a combined material 
degradation and geometrically nonlinear analysis is 
performed which is called a progressive failure analysis. 
In quantifying the failure state, the percentage of failed 
plies is used and defined as the number of plies with 
failures at each Gauss point divided by the product of 
the number of plies and the number of Gauss points. 

Rail-Shear Panel 

Rail-shear tests are frequently used to measure in- 
plane shear strength. The rail-shear specimen used in 
this analysis is the T300/976 graphite-epoxy cross-ply 
laminate reported on by Shahid 19 with both analytical 
and experimental data. 

The material properties of the T300/976 graphite- 
epoxy used in this analysis are based on the data 
available in Shahid. 19 These values include 
E u = 20.2 x 10 6 (psi), E 22 = 1.41 x 10 6 (psi), G n = 
0.81 x 10 6 (psi), and t> 12 = .29. Note that G 12 is a 
function of the stacking sequence. Based on Shahid, 19 
the material allowable values are 
X T = 220.0 x 10 3 (psi), X c = 231.0 x 10 3 (psi), Y T = 
6.46 x 10 3 (psi), Y c = 36.7 x 10 3 (psi), and in-plane 
shear strength T = 6.0 x 10 3 (psi). The 24-ply 
laminate has a stacking sequence of [0 6 / 90 6 ]^ . The 
nominal ply thickness is taken as 0.0052 inches. 

The specimen is 6-inches long and 1-inch wide. In 
order to represent the boundary conditions of the 
specimen in a rail-shear fixture, one edge of the 
specimen is firmly fixed, while on the other parallel 
edge, deformations are allowed parallel to the edge in the 
y-direction and restrained from motion in the re- 
direction. Upon loading, a displacement increment is 
applied along the latter edge. A finite element model of 
this specimen with 48 4-node elements along the length 
and 8 4-node elements along the width is shown in 
Figure 2. 

The structural response of the rail-shear specimen is 
studied. A nonlinear analysis is first performed to 
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understand the behavior of the rail-shear specimen 
without any material failures. Then progressive failure 
studies are performed on the model to evaluate the three 
failure criteria on predicting the failure of the panel. 
For the analyses, an initial displacement of 0.0005 
inches is applied to the left edge of the rail-shear 
specimen. The displacement increment for successive 
progressive failure analysis load steps is also 0.0005 
inches. 

A summary of the failure loads (first-ply failure and 
final failure) and the dominant failure mode type using 
the three failure criteria is presented in Table 1 . Load- 
deflection curves presented in Figure 3 compare all 
criteria with the experimental and analytical results of 
Shahid. 19 All progressive failure results agree very well 
with the experimental results. Figures 4 and 5 show 
the structural response at selected loads along the load- 
deflection curve using Christensen’s criterion. The left 
figures (Figures 4a and 5a) indicate the percentage of 
failed plies within an element, and the right figures 
(Figures 4b and 5b) depict the N stress resultant 

distribution. Dark regions indicate high values. These 
figures show the development of damage and its 
propagation to failure. The progressive failure analyses 
using other criteria show similar structural responses. 
The laminate fails along the edges of the panel where 
the shear loading is applied. 

Compression-Loaded Composite Panel 

The next example is a composite rectangular panel 
loaded in axial compression. The panel length is 20.0 
inches and the width is 6.75 inches. This panel is 
denoted Panel C4 in the experimental results reported by 
Starnes and Rouse. 20 

The panel is fabricated from unidirectional Thomel 
300 graphite-fiber tapes preimpregnated with 450K cure 
Narmco 5208 thermosetting epoxy resin. The material 
properties of the specimen are based on the data 
available from the test. 20 The values include E u = 

19.0 x 10 6 (psi), E 22 = 1.89 x 10 6 (psi), G n = 
0.93 x 10 6 (psi), and v l2 = .38 . From Engelstad’s 
paper, 6 the material allowable values are 
X T = 200.0 x 10 3 (psi), X c = 165.0 x 10 3 (psi), 

Y t = 1 1.74 x 10 3 (psi), Y c = 27.41 x 10 3 (psi), and 
in-plane shear strength T = 10.0 x 10 3 (psi). The panel 
used is a 24-ply orthotropic layup with a 
[±45/ 0 2 /±45/0 2 /±45/0/90]^ stacking sequence. 

The thickness of each ply is 0.00535 inches. 


The finite element model of this panel shown in 
Figure 6 has 40 4-node elements along the length and 
14 4-node elements along the width. The loaded ends of 
the panel are clamped by fixtures, and the unloaded ends 
are simply supported by knife-edge supports to prevent 
the panel from buckling as a wide column. 

The structural response of the C4 panel is studied. 
A fringe plot of the first buckling mode from the linear 
stability analysis is shown in Figure 7a and compared 
to the moire-fringe plot from the experiment 20 shown in 
Figure 7b. The results indicate that the first buckling 
mode from the analysis and the experimentally observed 
buckling mode shape are in agreement with each other. 
The first buckling mode has two longitudinal half- 
waves with a buckling node line at panel midlength. 
An initial geometric imperfection is formed by using 
the first buckling mode shape normalized by its 
maximum component, then scaled by 5% of the panel 
thickness, and added to the nodal coordinates. The 
initial displacement applied to the panel is 0.001 
inches. A displacement increment of 0.0025 inches is 
used for the first 10 steps in the progressive failure 
analysis. Then, a smaller increment of 0.001 inches is 
chosen for the next 5 load steps so the analysis can pass 
the buckling load. In the next 25 load steps, a 
displacement increment of 0.0025 inches is used for the 
analysis. Finally, a displacement increment of 0.001 
inch is used until ultimate failure of the panel is 
predicted. This adjustment of the applied displacement 
is done to reduce errors potentially introduced by not re- 
establishing equilibrium after material degradation. 

Progressive failure analyses using Hashin’s and 
Christensen’s criteria were performed on the C4 panel. 
Similar progressive failure analyses have also been 
performed by Engelstad, Reddy, and Knight. 6 The 
progressive failure results for the C4 panel are presented 
in Figures 8 and 9 using Hashin’s and Christensen’s 
criteria. The experimental and analytical out-of-plane 
deflections near a point of maximum deflection are 
shown in Figure 8 as a function of the applied load and 
correlate reasonably well. A comparison of analytical 
and experimental end shortening results as a function of 
the applied load is shown in Figure 9. The progressive 
failure results for end shortening also agree reasonably 
well with the experimental results. 

At some point in the progressive failure analysis, a 
dramatic change in the slope of the end shortening curve 
indicates an inability for the panel to support any 
additional load. This load level is designated as the 
analytical failure load, and the final experimental data 
point is called the test failure load. The final failure 
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loads predicted by both criteria are very close to each 
other. However, the first-ply failure (FPF) load of 
Christensen’s criterion is much lower than Hashin’s 
criterion. Both failure criteria predicted higher failure 
loads than the results from the experiment. Table 2 
provides a summary of the failure loads (first-ply failure 
and final failure) and the dominant failure mode type for 
both failure criteria and the test results. The dominant 
failure mode of Christensen’s criterion is fiber/matrix 
interaction, and the dominant mode of Hashin’s criterion 
is failure in matrix tension. Starnes and Rouse 20 
reported that the experimental failure mode was due to 
transverse shear effects at the nodal line of the buckle 
pattern near panel midlength as a result of coupling of 
large out-of-plane deflections and high transverse shear 
strains near the edges of the panel. Engelstad et al. 6 
analytically confirmed this behavior. The current 
progressive failure analysis capability does not include a 
transverse shear failure mode, thus, these analytical 
results cannot predict this failure mode, but still agree 
reasonably well with the test results. In the future, 
interlaminar stress predictions will be included in the 
COMET formulation to predict interlaminar failures. 

The structural response of the C4 panel at selected 
loads (see Figure 9) is given in Figures 10 and 1 1 using 
Hashin’s criterion. Figures 12 and 13 show the 
structural response at selected loads using Christensen’s 
criterion. Each set of figures displays the percentage of 
failed plies within an element, the N x stress 
distribution, and the out-of-plane deflection of selected 
steps in the progressive failure analysis. This sequence 
of figures shows the development of damage and its 
propagation to final failure. The results show that the 
damage is concentrated along the nodal line just as the 
test results indicated in reference 20. The out-of-plane 
deflection pattern exhibits large deflections (nearly three 
times the panel thickness at the buckle crest) and high 
local gradients. The high in-plane membrane stress 
resultants and their gradients near the buckle nodal line 
contribute to the failure propagation and final failure. 
The in-plane normal and shear stresses interact for the 
Hashin and Christensen criteria. Examination of the 
twisting moment resultant M shown in Figure 14b 

indicates that near the failure locations high M values 
exist. High M ^ values produce high shear stresses 
in the region where the failures occurred. This 
effect is compounded by high compressive stresses N x 
shown in Figure lib and contributes to understanding 
this type of failure. 


Concluding Remarks 

A progressive failure methodology has been 
developed, extended, and successfully implemented in 
COMET. The current progressive failure methodology 
uses C 1 plate and shell elements based on classical 
lamination theory. The progressive failure 
methodology accommodates various formulations in 
predicting failure such as the maximum strain criterion, 
Hashin’s criterion, and Christensen’s criterion. These 
different formulation methods are compared and assessed 
by performing analyses on laminated composite 
structures. Results indicate good comparisons with 
existing test data on two laminated composite panels. 

To accomplish a progressive failure analysis tool 
capable of predicting all types of failures, the current 
progressive failure methodology should be enhanced. 
The first enhancement is to extend the current 
progressive failure model such that it would include 
failure criteria to predict failure mechanisms involving 
interlaminar stresses. One option to accomplish this 
using the current progressive failure model is to 
calculate the interlaminar stresses (<7^ p yz O^) by 

integration of the 3-D equilibrium equations. Another 
option is to incorporate the new shell element by 
Tessler 21 which includes the interlaminar stresses as part 
of the element formulation. Once the interlaminar 
stresses are known, other failure mechanisms to predict 
debonding or delaminations could be added. Another 
enhancement is to extend the progressive failure 
analysis model so that static equilibrium is re- 
established after material properties have been degraded. 
This is accomplished by repeating the nonlinear 
analysis at the current load step until a converged 
solution exists. Such a capability permits the use of 
arbitrary step sizes during the nonlinear analysis and 
provides for an automatic step size control rather than 
fixed step size. 
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Table 1 Rail-Shear Panel: Summary of Progressive Failure Results. 


Failure Criterion 

First-Ply Failure Load 
(lbs.) 

Final Failure 
Load (lbs.) 

Dominant Failure 
Mode Type 

Christensen 

2105 

3964 

Fiber-Matrix Interaction 

Hashin 

2406 

4016 

Matrix Tension 

Maximum Strain 

2707 

4065 

Shear 

Test Results 19 

Unavailable 

4010 

Shear 


Table 2 C4 Panel: Summary of Progressive Failure Results. 


Failure Criterion 

First Ply Failure Load 
(lbs.) 

Final Failure 
Load (lbs.) 

Dominant Failure 
Mode Type 

Christensen 

18615 

22454 

Fiber-Matrix Interaction 

Hashin 

21778 

23526 

Matrix Tension/ 
Fiber Compression 

Engelstad et al. 6 
Maximum Stress 

22724 

25049 

Transverse Shear 

Engelstad et al. 6 
Tsai-Wu 

20927 

23464 

Transverse Shear 

Test Results 20 

Unavailable 

21910 

Transverse Shear 
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Figure 1 Typical Progressive Failure Analysis Methodology for ith step. 
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Shear Load, P (lbs.) 




6 in. 




Figure 2 Geometry of Rail-Shear Panel with 48 x 8 Finite Element Mesh. 



Applied Displacement, u (in.) 

Figure 3 Rail Shear Panel: Load vs. Deflection Response. 


10 




Max:= .65192E+03 



.58673E+03 
.52154E+03 
.45634E+03 
.391 15E+03 



.32596E+03 
.26Q77E+03 
.19558E+03 
. 1 3038E+03 
.65192E+02 


Min;= .OOOOOE+OO 


b) N xy Stress Distribution 


Figure 4 Rail-Shear Problem: Structural Response at First-Ply Failure (2105 lbs.). 
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a) Percentage of Failed Plies 


Max := 100% 



Min := 0 % 



b) N xy Stress Distribution 


Max:= .65192E+03 

.58673E+03 
.52154E+03 
.45634E+03 
.39115E+03 
.32596E+03 
.26077E+03 
.19558E+03 
.13038E+03 
.65192E+02 

Min:= .OOOOOE+OO 



Figure 5 Rail-Shear Problem: Structural Response at Final Failurer(3964 lbs.). 
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Out-of-Plane Deflection, w (in.) 

Figure 8 C4 Panel: Out-of-Plane Deflection Comparison. 



End Shortening, u (in.) 


Figure 9 C4 Panel: End-Shortening Comparison 
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a) Percentage of Failed Plies 


b) N x Stress Distribution 


c) Out-of-Plane Deflection 


Figure 10 C4 Panel: Structural Response at First-Ply Failure (21778 lbs.) - Hashin’s Criterion. 
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a) Percentage of Failed Plies 


b) N x Stress Distribution 


c) Out-of-Plane Deflection 


Figure 1 1 C4 Panel: Structural Response at Final Failure (23526 lbs.) - Hashin’s Criterion. 
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a) Percentage of Failed Plies 


b) N x Stress Distribution 


c) Out-of-Plane Deflection 


Figure 12 C4 Panel: Structural Response at First-Ply Failure (18615 lbs.) - Christensen’s Criterion. 



a) Percentage of Failed Plies b) N x Stress Distribution 


c) Out-of-Plane Deflection 


Figure 13 C4 Panel: Structural Response at Final Failure (22454 lbs.) - Christensen’s Criterion. 
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Figure 14 C4 Panel: Structural Response at Final Failure (23526 lbs) - Hashin’s Criterion. 
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